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DETERMINATION OF THRESHOLD FAILURE LEVELS OF SEMICONDUCTOR 
DIODES AND TRANSISTORS DUE TO PULSE VOLTAGES 

ABSTRACT 

Theoretical predictions of circuit failure in an Electromagnetic Pulse 
(EMP) environment require a knowledge of failure levels for each component of 
the circuit due to surge voltages or currents. For most circuits, the semi¬ 
conductor devices are the weakest elements with respect to such failure. This 
section presents the results of an extensive experimental program to determine 
pulse power failure levels of semiconductor junctions. Approximately 110 

different types of silicon diodes and transistors were studied with variations 
_li _i 2 

in junction areas from 10 to 10 cm and with widely varying junction geometries. 
Power levels of up to two kilowatts, with time durations of 0.1 to 20 microseconds, 
were applied to semiconductor junctions in both forward and reverse polarity modes. 
A semi-empirical formula, based on experimental data and on a simple thermal 
. failure model is given. From the formula one can make order-of-magnitude 

estimates of the failure level as a function of pulse length for many silicon 
diodes or transistors whose junction area is known. 

Introduction 


There have been extensive studies of nuclear radiation effects on 
electronic circuitry and systems for a variety of radiation environments. 

An area of increasing interest in the past few years is the effect of 
Electromagnetic Pulse (EMP) radiation. There are many phases in studies 
of EMP effects on electronic circuits. A major phase of such a study is 
the determination of the free-field EMP environment, and the transfer function 
between this environment and the system configuration under consideration. 
These results are used to predict the electrical signal impressed on the 
circuitry. Part’s of these studies have been published [l]-[3] and further 
consideration of this phase of the overall problem is not given in this 
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paper. A second major phase, the subject discussed in this section, is the 
study of effects once the induced circuit transients are known. Using the 
calculated amplitude and waveshape of the induced transient as an input 
parameter, one can calculate the circuit behavior by use of network analysis 
programs designed for computer use. As a principal part of the overall 
prediction of EMP effects on circuitry one would like to be able to make a 
reasonable estimate of whether or not the circuit will properly respond once 
the transient conditions have passed. Such predictions require a knowledge 
of the levels at which the circuit components will fail due to electrical 
stress. In general, the semiconductor components are the most susceptible 
to damage by transient voltages. Good composition resistors can withstand 
pulse powers of better than 10,000 times their power rating for microsecond 
pulses. Also voltage pulses high enough to cause semiconductor damage will 
not damage many resistors because the current is limited to very low values 
by the resistance of the device. 

There has been extensive work done on damage and failure mechanisms 
in semiconductors [U]~[9]. However, there is no general failure level 
information available for the extensive number of diodes and transistors 
used in present-day circuitry. Specifically, there is a particular lack 
of failure levels for microsecond and sub-microsecond pulses. Shielding 
circuitry in an EMP environment would become a formidable problem if the 
manufacturer's maximum ratings were used to establish a failure level for 
a device. A more logical solution is to determine the actual failure levels 
and shield accordingly. Also, if failure levels of components are known, then 
the problem of an EMP environment can be considered during circuit design and 
proper components chosen. The determination of failure levels is not only 
applicable to EMP problems but is also applicable whenever high transient 
voltages appear in a circuit whether the pulse origin is EMP, gamma induced, 
or a transient from within the system itself. 

The failure level program to be described had three main objectives: 

(l) to determine approximate failure levels for some 110 specific semi- 


1.2 






I 


I 

\ L 

; conductor diodes and transistors, (2) to determine whether one could obtain 

; a model or an empirical formulation which would enable one to calculate 

approximate failure levels for other devices, and (3) to calculate failure 
; levels for the devices tested at various pulse durations for which testing 

was not specifically accomplished. Such empirical formulations appeared to 
i be quite feasible based on earlier results of Davies and Gentry [10] for 

longer pulse durations. However, semi-empirical formulations of solutions 
to such a problem must be used with caution. Because of the wide variety 
and characteristics of semiconductor devices the analysis presented may not 
be valid in all cases. 


Testing Philosophy 

Failure levels for various pulse durations were desired for a large 
number of different types of devices. A large enough program to provide 
statistically valid data for each test was determined to be prohibitive, 

( Therefore, the testing program was based on small samples (generally 5 

to 10) of each device for each test in order to achieve a value for the 
general region in which threshold failure occurs. The confidence in the 
failure levels so obtained is greatly increased by providing a general 
theoretical framework in which to evaluate the experimental results. A 
thorough inter-comparison of the different devices is made based on their 
common modes of failure. Since there is an infinite number of conditions 
under which a transistor may be operating (i.e., different circuit connections), 
the transistor failure was found for each terminal pair with the third terminal 
open. Predominantly the emitter-base junction was tested. No steady state 
power or bias was applied to the transistor during tests. Failure analysis 
can then be correlated with circuit code predictions for the power level 
between each specific terminal pair. 

If a device is tested in an unbiased condition, the behavior of certain 
microscopic processes will change. A number of junction phenomena, such as 
second breakdown [ll]~[lU], depend on the total combination of transistor 
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conditions. The present computer circuit code models are also not of 
such a nature to take into account the many complex physical phenomena 
occuring. Therefore, the failure level of some devices under certain 
operating conditions might well he significantly lower than the following 
analysis predicts. There has not been-sufficient comparison between 
failure levels obtained in the manner described to those obtained in 
specific circuit configurations to be able to state that such analysis is 
always valid. However, it is not expected that the failure levels will be 
orders of magnitude higher because the analysis is based on the formation 
of microplasmas and hot spots within the junction. The presence of local¬ 
ized molten areas within the junction is nearly certain to cause permanent 
changes in the junction. Only in the cases of very small area microplasmas 
has full recovery been reported. 

Junction Failure Modes 


Because the present study involvement is concerned with high amplitude 
voltage pulses between two terminals of an inactive device, it is possible 
to examine the failure mode of a single p-n junction whether or not the actual 
device is a diode or a multiple junction device. However, for a more comprehensive 
study of possible failure modes, it is necessary to be aware of all processes 
that may be peculiar to the device under consideration. Such processes can depend 
on circuit conditions, the number and types of p-n junctions, and on general 
geometrical and construction details. 

The failure of a solid state junction device is considered to have 
occurred if the parameters of any p-n junction have been seriously degraded 
or if the device junction has become an open or a short circuit. In a more 
comprehensive study one might wish to study the degret of the degradation 
of parameters of the device such as gain, junction capacitance, reverse 
saturation current, breakdown voltage levels, etc., below a certain pre- 
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determined level. Although the actual failure of the device is due to an 
opened or shorted junction, it is convenient for analysis to discuss the 
failure in terms of the mechanisms which eventually cause the open or short 
to occur. It is also convenient to divide the discussion between effects 
due to the application of forward or reverse voltages. 

The principal breakdown mechanisms for reverse voltages on a single 
p-n junction are: 

(1) surface breakdown around the junction, and 

(2) internal breakdown through the junction within the body of 
the device. 

The problem of surface breakdown over the junction has received 
considerable design consideration [10] in more recent years. With the 
understanding of theoretical models based on the solution of Poisson's 
equation with different geometrical boundary conditions, junctions have 
been successfully designed and built which exhibit body breakdown prior to 
any surface breakdown phenomenon. The problem of theoretically predicting 
surface breakdown is a very difficult one since it depends upon so many 
parameters such as geometrical design, doping levels near the surface, 
lattice discontinuities on the surface, and general surface conditions. 

The theoretical models which have been used are generally limited to the 
calculations of field-gradients at the surface under homogeneous crystal 
conditions and approximate geometrical boundaries which lend themselves to 
theoretical calculation. Under pulse conditions, the problem is even more 
formidable. Davies and Gentry state, "Unfortunately the transient energy 
which can be dissipated during surface breakdown is both unpredictable and 
appreciably lower than that which can be absorbed within the body of the p-n 
junction device." 

The actual destruct mechanism of such a surface breakdown is usually to 
establish a high leakage path around the junction, thus nullifying the junction 
action. The actual internal junction itself is not necessarily destroyed as 
ha3 been shown since re-etching the surface can return the diode to normal 
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action. In internal tody "breakdown the destruct mechanism apparently 
results from changes in the junction parameters due to the high temper¬ 
atures locally within the junction area. These temperatures can "be of 
such magnitude that alloying, or diffusion of the impurity atoms, occurs 
to such an extent that the junction is either totally destroyed or its 
properties drastically changed. Figure 1.1 shows four mechanisms for the 
movement of impurity atoms through a lattice. 
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(d) By means of an 


Figure 1.1 Mechanisms of Movement of Impurity 
Atoms through a Lattice. 


Of the four mechanisms depicted, (a) direct interchange, (b) cooperative 
interchange, (c) interstitial, and (a) adjacent vacancy, the last two are 
predominant in silicon. 

High temperatures and strong applied fields as well as the density gradient 
aid the movement through the lattices. 




In many cases both for forward and reverse voltages the current is 
sufficiently high and localized to cause melting at hot spots [4]-[9] 
within the junction. Such action can result in a resistive path across 
the junction which masks any other junction action. 

For long term heating at temperatures well below the melting level 
impurity diffusion can be an important mechanism contributing to device 
degradation. For very short pulses hot spot formation is more likely the 
predominant cause of failure. 

Junction failure can also occur due to dielectric breakdown. The large 
avalanche current can form a path for an arc discharge to occur. This can 
result in a puncture through the junction with an actual pin hole being formed, 
Usually a junction short results. This effect will be described as a puncture 
rather than "punch-through" reserving the latter term for a depletion width 
phenomena occurring in multiple junction devices. 

For forward voltages surface breakdown is not expected to be a problem. 

The field gradients on the surface across the junction are normally quite 
small since the junction is in a highly conducting state. 

The failure of the junction for forward voltages is then expected to 
be primarily due to the temperature rise due to the high currents passed 
during the pulse. The destruct mechanism being mainly the change in junction 
parameters due to the high temperatures as discussed for reverse voltages. 

To illustrate a condition in which the operating conditions are important, 
the .breakdown mechanism occurring in transistors called punch-through will 
briefly be discussed. If the base region is of higher resistivity than the 
collector region, the depletion region of the collector-base junction extends 
mainly into the base region. As the collector voltage is increased, the 
depletion region extends further and further into thv base region. If the ■ 
base region is very thin, then the collector voltage may be high enough to 
cause the depletion region to extend entirely through the base and make contact 
with the emitter-base depletion region thus effectively causing a short across 
the base region. This will occur if the collector voltage necessary to extend 
the depletion region through the base is smaller than the voltage required to 
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initiate avalanche breakdown. Again, if the current is limited there is no 
permanent damage done to the transistor due to the punch-through. When the 
collector voltage is decreased the transistor will return to normal action. 

The pulse rise time also may affect the breakdown characteristics of 
a junction. Experiments by Agatsuma, Kohisa, and Sugiyama [15] on the collector 
to emitter breakdown with the base open showed that the energy required to 
initiate second breakdown was much less with fast rising pulses. One theory 
is that if the voltage rise is slow, the breakdown occurs over the whole 
junction area. If the voltage rise is ve*y fast, it appears that very 
localized breakdowns occur within the junction which then carry very large 
currents. Ioffe and Regel [l6] give the conductivity of molten silicon as 

4 -1 -1 

about 10 ohm - cm which is about 30 times greater than the solid at 
the melting point. Thus,• appreciable currents can be conducted by small 
areas of the junction. Localized hot spots occur and the temperature rise 
can be great enough to cause permanent damage. This can occur at pulse 
energy levels much below the energy levels required to cause damage with 
slow rising pulses. 

Considering nanosecond rise-time pulses, it could be expected that the 
reverse voltage breakdown level would become appreciably higher since the rate 
of rise could be appreciably faster than the time required to initiate breakdown. 
Experimental results published by Portnoy and Gamble [17] for collector- 
emitter breakdown with the base open show that voltages above breakdown 
can be applied for fairly long times (tens of microseconds) without break¬ 
down. If the pulse continues for longer times, the junction breaks down 
and the voltage drops to the calculated breakdown voltage. 

Also, depending on the pulse rise time and width, other parameters 
such as the skin effect, non-homogenous local fields, emitter crowding 
and ion migration may affect the failure levels. 


Modeling for Junction Failures 

Although many different microscopic mechanisms occur, it has been 
found that most of these failure mechanisms are linked primarily to the 
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junction temperature. Therefore, the theoretical treatment of the problem 
can be reduced to a thermal analysis. The pulse times considered in this 
study are 0.1 to 20 microseconds. For these short times, the boundaries 
of the bulk material and the thermal heat sinks do not have the effect they 
would at longer times. The junction temperature can be adequately approximated 
from simple linear heat flow theory. Although more exact thermal models have 
been devised [10], [l8j-[20] to account for specialized boundary conditions, 
the attempt here is to use an intermediate model that is reasonable for a 
wide variety of devices. In this manner the experimental results for dif¬ 
ferent devices and pulse duration can be readily compared. There is no doubt 
that a better model could be made for any one specific device if sufficient 
information on the junction parameters such as geometry, thickness, doping, 
resistivity, and the geometry of the device as a whole were available. Inves¬ 
tigations were not made during this study to attempt to relate failure to any 
of these other parameters. 

Microplasmas and hotspot information are presently being investigated and 
results reported by many authors. As better three-dimensional models of hot¬ 
spot formation are made, better estimates of junction burnout as a function 
of pulse width may be made. One of the principal problems of three dimensional 
solutions is the need for more exact information about the geometry of the 
devices. Generally such information is not readily available. Other models 

of the temperature rise of hotspots may give a slightly different time 

- 1/2 

dependance. It is felt that the time dependance of t derived here is 

appropriate and sufficiently accurate for the pulse times considered. The 
two dimensional model derived by Takagi and Mano [20] gives a maximum current 

-lA -1/2 

time dependance of t and hence a power dependance of t . The worst 
case as far as achieving high temperatures in the junction is when one 
considers that all of the power dissipated in the device occurs in the junction. 
This corresponds physically to the situation where a high-voltage pulse of 
reverse polarity is applied to a junction with a high reverse voltage break¬ 
down. When the avalanche breakdown occurs, almost all of the applied voltage 
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is dropped across the junction, ana~Or*~y.a sir xl percentage is dropped across 
the bulk material. For such a thermal model, the junction failure conditions 
(derived in Appendix A) are 


P/A =‘\j TrKpC [T - T.] t~ l/2 , 
’ p m 1 


( 1 ) 


where P is the power; A, the junction area; k , the thermal conductivity; p, the 
density; C^, th 2 specific heat; T ffl , the failure temperature, T^, the initial 
temperature; and t, the time. This equation is of similar form to that suggested 
by Davies and Gentry [10]; its simple form leads to many distinct advantages in 
general use. Equation (l) is plotted for three specific cases in Fig. 1.2: 

Case I - Heating from approximately room temperature of 25°C to 
675°C which is the model failure temperature suggested 
by Davies and Gentry and is also within the temperature 
range discussed by other authors [4], [9]. 

Case II - Heating from 25°C to the approximate melting temperature 
of silicon of l4l5°C. 

Case III - Identical to Case II except that the currents are assumed 
to pass through localized hot spots which account for 
one-tenth the total junction area. This choice provides 
a convenient reference when intercomparing data since the 
curves for Case II and III are exactly separated by one order 
of magnitude. 

The form of the solution has many advantages such as: (l) It is a straight 
line on log-log plots ar.d hence allows a quick estimate of failure levels by 
plotting the power versus pulse time duration for a known junction area. This 
is in contrast to more exacting models whose solution is usually of the form 
of an infinite series and require more exacting information about the specific 
geometry of the device. (2) The only physical parameter needed for a specific 
device is the junction area. (3) It is simple to use as a semi-empirical 
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equation to match specific experimental data by keeping the appropriate 
slope of -1/2 and adjusting the constant to best fit the data. It is 
possible to determine the constant from the curves shown by merely reading 
off the value of the power per unit area at the one microsecond pulse duration 
provided the curves are plotted in units of kilowatts/sq.cm, and microseconds. 

(1+) It allows an intercomparison of devices of different junction areas and 
various junction geometries by plotting the power per unit area for failure 
for the different devices all on a single plot. 

In order to utilize these advantages, the solution was also used to 
compare data in those cases for which it was not used as an appropriate model. 

Two examples of this are (.1) devices where the reverse breakdown voltage is 
small and (2) when the greater percentage of the applied forward voltage is 
dropped across the bulk material rather than the junction. In such cases 
the actual temperature of the junction is lower than the model implies. 

However, for the pulse durations under consideration, the calculated slope 
is still approximately valid and when comparison is made between the theoretical 
curves and the experimental data, the experimental curves will prove to be 
higher. If the reverse failure levels are used for all cases it merely gives 
a slightly conservative estimate for the forward voltage mode. At very short 
pulse durations (sub-microsecond) extremely large powers are required to obtain 
failure. Hence in many cases large currents are required with the result that 
significant amounts of power are dissipated throughout the bulk semiconductor. 

In such short times the heat transfer to the junction from points away from the 
junction is very small. For such cases the temperature of the junction approaches 
a dependenc. more like l/t rather than l/>/ir. However, within the time scales 
studied the deviation was not sufficient to justify a separate model. Again, . 
the failure levels presented are on the conservative side for very short times. 

An inherent limitation of the power failure model is the fact that actual 
failure modes are not one dimensional heat flow processes, and the simple 
theory cannot be expected to give exacting answers. However, certain of the 
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assumptions tend to have a cancelling effect. For example, phase change 
energy considerations and bulk heating (for devices where the reverse break¬ 
down voltage is small or in the case of forward applied voltages) would 
raise the power curve of Fig. 1.2, while hot spot formation over less than 
one-tenth the total junction area would lower it. The experimental data 
supports the use of the theoretical estimates even with such a wide range 
of junction geometries and areas, failure mechanisms, and thermal boundary 
conditions. 

To perform a complete theoretical failure analysis it is required to 
know the power being delivered to the semiconductor junction. This inv.olves 
calculating the voltage and current levels for the junction. For this 
program, the power levels for failure were determined experimentally. 
Therefore, a thorough discussion of models for junction voltage and current 
levels will not be made. A very brief discussion of some of the basic 
approaches to voltage and current calculations is given in Appendix A. 

Experimental Conditions and Results 

All of the devices were tested to determine the threshold failure level 
of semiconductor junctions due to a high electrical stress. The general 
procedure consisted of step voltage pulsing the junction until failure or 
degradation of the device occurred. In the majority of the failures, a 
single pulse reduced the DC Beta and the zener voltage to zero with the 
device assuming a resistor characteristic. Pulse power levels of up to 2 
kilowatts (a few peaks to 8 kilowatts) were applied with pulse times of 
100 nanoseconds to 20 microseconds. The risetime of the applied pulse was 
approximately 20 nanoseconds. For each pulse, the waveform of the current 
through and the voltage across the device were recorded photographically. 

The photographic data was reduced to show the power level and the pulse 
duration applied to the terminals of the device. Figures 1.3 and 1.U show 
typical photographs for a reverse polarity pulse applied to a junction. 

The pulse in Fig. 1.3 caused no damage. The pulse in Fig. 1.1* caused the •' 
junction to fail during the time that the pulse was applied. The applied. 
















pulse time was approximately the same in hoth figures. Figure l.-'t shows that 
when the junction failed the voltage across the terminals dropped by 
approximately the value of the reverse breakdown voltage of the junction. 

An attempt to determine the mode of failure on each device was not 
made. However, a number of devices were examined after failure. In each *. 
case where the damage was evident, it appeared that localized melting had 
occurred. In some devices sufficient melting had taken place that a flow 
of melted silicon could be observed. From these results and the information 
of other investigators £1*]—[9]» it is believed that the predominant failure 
mode, for these microsecond and sub-microsecond pulses, is localized melting 
across the junction. The melted regions form resistive paths across the 
junction which shorts out or masks any other junction action. 

Since most devices tested required greater power for failure in the 
forward voltage mode, the experimental measurements were concentrated primarily 
on the reverse voltage mode. 

Figures 1.5 and 1.6 are power versus pulse duration plots for two specific 
diodes. The line on each plot is a line with a slope of -1/2, which best 
fits the data taking into account both failure and no failure points and 
uncertainty in the data for each point. The purpose of the figures is to 
give an indication of the actual power levels and the spread of points which 
caused failure. The plots illustrate that the time dependance is based on 
the total semi-empirical analysis rather than from the few data points obtained 
for each device. 

It was impossible to achieve a variation in burn-out time for several 
diodes tested for the reverse polarity voltage case. Normally a variation 
in burn-out time was achieved by reducing the amplitude of the applied 
voltage pulse. Also, when many of these diodes failed, the power level was 
below that previously dissipated without damage. A possible explanation for 
these results is as follows: Because of the geometrical uniformity of the 
diode, normally a large area of the junction is conducting and carrying the 
reverse current. Therefore, larger powers can be dissipated without failure 








Fig. 1.6. Experimental data points for failure of the 
anode-cathode Junction of lNh59 diode for reverse polarity 
voltage pulses. 







as local hot spots are not as likely to form. However, as the temperature 
of the junction increases and a certain voltage level is reached, small 
differences in the junction cause sufficiently high local field gradients 
that a dielectric breakdown occurs or the current is channeled into localized 
areas. When this happens, a local hot spot quickly develops and the device 
burns out. If this type of mechanism is occurring, longer burn-out times at 
smaller power levels could be achieved by using a pulser which contained 
current limiting circuitry. In a few of these devices the oscilloscope 
photographs would indicate that junction failure had occurred; however, 
when tested, junction action would still be present although the device 
would be very leaky. Possibly very microscopic breakdowns form a high 
resistance (low compared to the normal reverse resistance) path instead 
of a complete short across the junction. 

Figures 1.7 to 1.10 are power versus pulse duration plots for four transistors. 
Figures 1.7 and 1.8 show both forward and reverse voltage failure points for the 
base-emitter junction. Figure 1.8 also includes damage levels for the collector- 
emitter terminals. It should be noted that while the collector-emitter terminals 
involve txo junctions, the failure level lies within the general failure region 
for a single junction. For the collector-emitter case it was found that either 
one or both junctions failed. 

Figure 1.11 shows a composite of'the failure points for six types of 
semiconductor diodes inter-compared by plotting the power per unit area for 
failure versus pulse duration. Figure 1.12 shows a^omposite of base-emitter 
junction fr.'lure points for eight types of transistors inter-compared on a 
power per unit area basis. Figure 1.13 shows a composite of ten large area 
diodes for which failure points were not achieved. It was impossible to cause 
failure in some of the large area diodes because of the current limitation 
of the pulser. Because of the small reverse voltage breakdown for many of 
these diodes, large currents of the order of 30 to 1*0 amps can be carried 
without damage. Several other diodes could not be damaged because of the 
voltage limitations of the pulser (1KV maximum) and the high peak inverse 
voltages of the devices. The reverse voltage breakdown levels of these diodes 
were either not reached or were only slightly exceeded with the microsecond 




Fig. 1.7. Experimental data points for failure of the 
base-omitter Junction of a 2N2222 transistor for forward 
and reverse polarity voltage pulses. 



Fig. 1.8. Experimental data points for failure of the 
base-emitter and collector-emitter Junctions of a 2N1132 
transistor for forward and reverse polarity voltage pulse 
















fig. 1.11. A composite of experimental data points for the 
failure of the anode-cathode junction for six diodes. 
Inter-comparison is made by plotting power per unit area 
versus pulse duration. Also shown are the theoretical 
failure curves. 



Fig. 1.12. A composite of experimental data points for the 
failure of the base-emitter junction for eight transistors. 
Intercomparison is made by plotting power per unit area 
versus pulse duration. Also shown are the theoretical 
failure curves. 











pulses. Therefore, no damage to the diodes was possible with the low 
currents that resulted. In Fig. 1.13, the power per unit area levels 
shown are generally all below the power per unit area levels for which 
failure was observed in other similar devices. Th^re is no certainty 
that these devices would not be damaged if these higher power levels had 
been reached. 

Some of the plots of power per unit area versus pulse duration 
may be somewhat in error because of the uncertainty in determining the 
junction area. Approximately one-third of the area values for the tested 
devices were obtained from the manufacturer; the remaining were calculated 
from physical measurements of the devices. It is felt that no serious 
misinterpretation could result in an errant junction area calculation since 
a departure from an actual area value by a factor of two would not result 
in a change any greater than the spread of points already existing. 

Figure l.lU shows a composite of the failure points for a 2N2222 
transistor tested under this program and data from a paper by Davies and 
Gentry [10] who gave failure points for some longer duration pulse lengths. 
The data is in general agreement for nearly six orders of magnitude in 
pulse duration. 

Each of the damage lines similar to those shown in Figs. 1.5 to 1.10 
were converted to a power per unit area basis. The constant for these lines 
was then averaged to find the best fit for all of the devices tested. The 
resulting equation is 

P/A = 510t~ 1/2 (2) 

where, for the time given in microseconds, P/A is in units of kilowatts/cm 2 . 
Shown in Figure 1.15 are lines obtained by averaging the diode and transistor 
results separately. The equations for these lines are: 

P/A = 55to _1/2 (3) 

and 

P/A = Vr0t- 1/2 (1+) 
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respectively. By using the above equations and the approximate junction area, 
a very good order of magnitude estimate of the failure level for most silicon 
diodes and transistors can be calculated. 

Table 1.1 gives values, for nineteen specific diodes, of the pulse 
power for a one microsecond pulse required to cause junction failure. 

The values range from 73 watts for a 1N2929A, a small silicon tunnel 
diode, up to 2.1 kilowatts for the 1N711A, a 7 volt zener diode. 

Table 1.2 gives similar values for eighteen specific transistors. The 
values range from 45 watts for a 2N930, a ..small planar KPN transistor, up 
to 2.3 kilowatts for an R22705638, a large area NPH transistor. 

Table 1.3 gives some typical experimental values of terminal voltage 
and current that caused failure for several diodes. The values range from 
high voltages of 1000 volts and small currents- of 1 amp for high reverse 
breakdown diodes to low voltages of 22 volts and high currents of 36 amps 
for a low voltage (6.2 volts) zener such as the 1N753A. Also shown are 
some typical voltage and current values for forward polarity tests. For 
most of the forward tests high ciUTents of up to UO amps could be drawn 
without failure for large area diodes. 

Table 1.4 shows similar values of voltages and currents for several 
transistors. 


Application of Results 

Computer circuit analysis codes such as CIRCUS, SCEPTER, NET-I, etc. 
[l8], are used to correlate the experimental failure analyses with specific 
circuit conditions. The elements of the circuit or system are arranged in 
nodal fashion to reproduce the circuit in the computer program. In the 
operation of these codes, equivalent circuit models for the semiconductor 
devices are available frc®. internal programming or as external inputs. 

The experime tal equivalent circuit parameters for the specific semiconductor 
devices and the calculated amplitude and waveshape of the transient signal 
at some node in the circuit complete the input information necessary for the 
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DIODE TYPE 


W6807WS 

W6807ZV 

1N711A 

1N823 

G129 

1N816 

1N981B 

IN753A 

1N702A 

UTl»59A 

IW482A 

1N540 

IH1095 

1R253 

1N967B 

1N537 

DHD936 

MC357 

1N2929A 


TABLE 1.1 Pulse Power 


Reverse 


MOMINAL FAILURE POWER 


(WATTS) " ) |i SEC. PULSE 

2200 

22C0 

2100 

1800 

1700 

1500 

lltOO 

1200 

1000 

960 

960 

930 

880 

860 

730 

510 

ll*0 

71* 

73 

one Microsecond Pulse required for 
r Junction Failure. 
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TRANSISTOR TYPE 


NOMINAL FAILURE POWER 
(WATTS) „ i y SEC. PULSE 


H22705638 

2300 

R227075 1 »7 

1900 

MIS 17181/1-1 

1500 

2N657A 

1070 

2N1116A 

980 

2Nl»98 

800 

2N657 

620 

2N335 

550 

2N33t> 

550 

2N1132 

500 

2N1893 

1*00 

2N336A 

31*0 

2N 699 

250 

2Nl61*2 

130 

2N2222 

110 

2N736 

100 

2N927 

96 

2N930 

2*6 


TABLE 1.2 Pulse Power for one Microsecond Pulse Required for 
Reverse Polarity Base-emitter Junction Failure. 
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DIODE TYPE TYPICAL REVERSE FAILURE TYPICAL FORWARD LEVELS 



LEVELS (~0.5 

- 2u SEC. PULSE) 

(-0.5 - 2m SEC. PULSE) 

(No Failure Unless Noted) 


Terminal 

Current 

Terminal 

Current 


Voltage 

(Volts) 

(Amps) 

Voltage 

(Volts) 

(Amps) 

W6807WS 

1000 

2 

5 

39 

1U711A 

72 

28 

15 

36 

1N823 

29 

3 6 

25 

3k 

G129 

ho 

28 

h 

10 

1N816 

710 

2 

19 

39 

1N981B 

175 

8 

29 

30 

1N753A 

22 

36 

11 

38 

1H702A 

20 

30 

111 

33 

lNl*59A 

1000 

0.9 

15 

36 

lNl*82A 

800 

1 

16 

36 

1N1095 

800 

1..6 

9 

35 

1N967B 

50 

ll» 

12 

31 Failed 

1N537 

1000 

1 

32 

1*0 

DHD936 

120 

1.5 

10 

16 Failed 

MC357 

160 

0.5 

15 

6 Failed 

1N2929A 

10 

8.5 

7 

G Failed 


TABLE 1.3 Values of Terminal Reverse Voltage and Current That Caused 
Junction Failure. 



TRANSISTOR TYPE 

TYPICAL REVERSE FAILURE 

TYPICAL FORWARD 

LEVELS 


LEVELS (-0.5 - 

2pSEC. PULSE) 

(-0.5 - 2pSEC. PULSE) 

(No Failure Unless Noted) 


Terminal 

Current 

Terminal 

Current 


Voltage 
(Volts) 

(Amps) 

Voltage 

(Volts) 

(Amps) 

R22705638 

55 

35 

12 

36 

2N1U6A 

70 

14 

28 

30 

2N498 

55 

24 

20 

32 

2N657 

62 

14 

22 

39 

2N335 

70 

6 

50 

15 

2N1132 

80 

8 

16 

16 Failed 

2N1893 

24 

14 

15 

14 

2N699 

68 

1.6 

— 

— 

2N1642 

130 

1 


— 

2N2222 

38 

3 

10 

14 Failed 

2N736 

10 

6 

~ 

- 

2N927 

320 

0.25 

~ 

— 

2N930 

30 

1.5 

20 

19 Failed 

TABLE 1.4 

Values of Terminal Reverse Voltage 

Base-Emitter Junction Failure. 

and Current That 

Caused 



computer program to analyze the response of the circuit. The computer 
calculates the time history of current through and voltage across any two 
terminals of semiconductor devices and then plots a time history of the 
power dissipated. The calculated pulse length and power information form 
the circuit code is correlated with experimental semiconductor failure curves 
to estimate failure of the semiconductor device due to the specific transient 
"being studied. 

The circuit codes and failure levels are also useful in circuit design. 
Specific environments can "be considered during the design phase and proper 
components chosen. Additionally, circuit reliability can be increased by 
incorporating protective circuitry where weak components exist. 

Conclusions 


It is acknowledged that while the approximate theories presented in this 
discussion are able to correlate a large body of experimental data, the phenomena 
involved are of such a complex nature that it cannot be presupposed that such a 
theory is applicable in all cases. To be useful, order of magnitude estimates 
must be used judiciously within the framework for which they were obtained. 

The mechanisms involved in microplasma and hot spot formation can be significantly 
affected by steady state operating conditions, structural geometry and to some 
extent by transient pulse waveshape. At longer and shorter pulse lengths the 

model described is not as appropriate. Even in the sub-microsecond times presented 

- 1/2 

here there is evidence that the curves are not holding to a strict t dependence. 
This is partially due to the basic one-dimensional model and partially due to the 
large amounts of power being dissipated in the bulk material rather than at . 
the junction. Direct surface effects and breakdown are also not covered in 
this analysis. Surface effects do not necessarily require high junction 
temperatures to cause gain degradation [22]. Veloric and Prince [23] state: 
"Experiments show that devices which show surface breakdown will collapse 
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at power levels which are orders of magnitude "below that observed for devices 
in which body breakdown is observed." 

Application of the analysis presented can produce satisfactory results 
providing the exceptional cases are recognized. 

Summary 


A simple thermal model was used to obtain estimates of semiconductor 
junction failure and provides the time dependence for a semi-empirical 
equation. This equation represents a best fit to experimental results from 

over 1000 devices of approximately 110 different types. The semi-empirical 
- 1/2 

relation is P/A = 510t , where, for the time given in microseconds, P/A 

is in kilowatts/cm . Actual failures on all types of devices tested were 
within better than one order of magnitude to the curves predicted by the 
thermal model and even closer to the semi-empirical equation. The semi- 
empirical failure equation is useful in that it allows for quick estimates 
of junction failure both for circuit design and for studies of electrical 
transients on existing circuitry. The equation can be used on most all 
types of silicon diodes and transistors whose junction area is known. 
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SECTION 2 


r 


! 


ESTIMATES OF SEMICONDUCTOR FAILURE 
DUE TO MULTIPLE VOLTAGE PULSES 

Previously, we have made estimates of semiconductor failure for a 
single applied voltage pulse. These estimates were based on thermal 
failure using a simple model to determine the temperature of the semi¬ 
conductor junction as a function of the pulse power and duration. 
Semiconductor devices were also tested for failure in the laboratory 
using a single voltage pulse. However, under certain conditions, one is 
interested in the failure level from a voltage wave which contains 
multiple pulses. This can arise directly from the induced EMP wave which 
may have a waveform as shown in Fig. 2.1 or it may arise from the circuit 



response due to a step input. Such a circuit response may be of the form 
shown in Fig. 2.2, 
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Our problem, then, is to estimate how such waveforms change the 
failure levels previously predicted for a single pulse. First, w». will 
discuss some failure estimates based on our original failure model. 

Secondly, we will discuss some of the failure mechanisms which may be- 
occurring and how these might cause the actual failure level for multiple 
pulses to differ from the predicted levels. 

Since semiconductor devices are not bilateral devices, we need to 
consider several po.s.'.ible cases separately as representative of what the 
operating conditions might be. 

Operating modes: 

1. Diode with high reverse voltage breakdown—First pulse 
in forward direction—very small reverse conduction. 

2. Diode with high reverse voltage breakdown—First pulse 
in reverse direction—very small reverse conduction. 

3. Diode with low reverse voltage breakdown or Base-Emitter 
junction of transistor—First pulse in forward direction- 
appreciable reverse conduction. 

1*. Diode with low reverse voltage breakdown or Base-Emitter 
junction of transistor—First pulse in reverse direction- 
appreciable reverse conduction. 

Our thermal model is bilateral since it is based on the power dissipated 
in the junction; therefore, we need to know the power dissipated in the 
junction for a given applied voltage waveform. Figure 2.3 shows the general 
waveforms one could expect for the four cases given above. 

For cases 1 and 2, there is a very small amount of power dissipated 
during the time the device is under a reverse voltage pulse. Fbr cases 
3 and h, current is flowing all the time since we have assumed that the 
reverse breakdown voltage is small compared to the applied voltage. For 
most cases, more power is dissipated directly at the junction when the 
device is conducting in the reverse direction. This occurs because of the 
much larger voltage drop across the junction in the reverse direction. 
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APPLIED PULSE 
VOLTAGE 




POWER WAVEFORM 
ACROSS JUNCTION- 
Case 2 


POWER WAVEFORM 
ACROSS JUNCTION- 
Case3 




POWER WAVEFORM 
ACROSS JUNCTION- 



Case 4 


Figure 2.3 










To simulate the above power waveforms in our theoretical model 
11 calculate failure for the following types of cases as shown in 
5. 2.1*. 






DOUBLE PULSE WITH 
COOLING PERIOD 

'ALTOPULSE DURATION 


DOUBLE PULSE WITH 
COOLING PERIOD EQUAL 
TO TWICE PULSE DURATION 


Figure 2.4 






For all cases considered, we willi assume „an equal power level for the 
pulses and will plot the results in terms of the individual pulse duration 


SINGLE PULSE: 

The failure levels for a single pulse were 1 "'previously calculated 
from equation 1, 


P 

A 


4tcpC [T -T. ]t -1 ^ 2 


( 1 ) 


where P is the power; A, the junction area; k, the thermal conductivity; 
p, the density; C^, the specific heat; T^, the failure tempeVature; 

, the initial temperature; and t, the time. For the cases where the 
failure temperature was assumed to be the melting temperature and (a) 
melting occurs over the entire junction and (b) melting occurs in local 
hot spots of one-tenth the total junction area, equation 1 becomes 
respectively: 


| = 1810 t" 1/2 (2) 

and 

1 = 181 (3) , 

} 

These equations, along with the averaged curves f~om actual experi¬ 
mental single-pulse failure data on diodes and transistors, are shown 
in Fig. 2.5. 

MULTIPLE PULSE: 

For multiple pulses, the worst case (lowest failure level) occurs 
when there is essentially no time between pulses for cooling. This is the 
type of pulses we see in cases 3 and of Fig. 2.3. For analysis purposes, 
we will simulate such waveforms by multiple square wave pulses as shown in 
Fig. 2.k (b, c) for the double and triple pulse. The use of square pulses 
will give a slightly lower failure level than would actually occur for 
sinusodial waveforms. 


** These averages change slightly as more devices are tested. The analysis 
is still valid and t.he curves are only shifted slightly in a linear manner 
depending on the averaged value. 
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DOUBLE PULSE AND TRIPLE PULSE: (No cooling period) 

Multiple pulses with no cooling between pulses can easily be calculated 
since the failure level is just the single pulse failure level for a 
single pulse whose duration is the sum of the multiple pulse durations. 

We can derive a general expression for the failure level for any number 
N of multiple pulses. 

From equation 1, we have that the temperature is 


P t V2 

[T f T.] = -- . (U) 

Ak p c p 


For a single pulse of duration t we have. 


P t 1/2 



(5) 


Then for a series 
equal power levels, we 


of N multiple pulses each of duration t and of 
obtain 


(T -T. ] 



( 6 ) 


and that 



P _ 
A ” 

N" 1/2 

{v^c" (T -T. ] T~ 1/2 } 

p mi 

(7) 

Equations 

2 and 3 

for multiple pulses are then 




P 

A “ 

»- 1/2 1810 

(8) 

and 







P _ 
A 

r l/2 181 t -l/2 

(9) 
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The constant K for the equation 





( 


P 

A 


- 1/2 


( 10 ) 


for the theoretical estimated and for the averaged experimental data for 
transistors and diodes is given in Table I. The curves for single, double, 
and triple pulses for the above cases are given in Figs. 2.6 (Theoretical), 
7 (Diodes), and 8 (Transistors). 


TABLE I 

Number of 

Value of K of Eouation 10 ! 

Pulses 

Theoretical 

Curves 

Diodes 

Transistors 

1 

1810 

181 

560 

310 

2 

1280 

128 

396 

219 

3 

10U5 

105 

323 

179 

h 

905 

91 

280 

155 

5 

809 

81 

250 

139 

6 

739 

7l* 

229 

127 


Estimated failure levels for double and triple pulses for a 2N2222 
transistor are shown in Fig. 2.9- These estimates are based on the actual 
experimental failure level for single pulses and on the above analysis. 
DOUBLE PULSE: (with cooling period) 

In order to make failure estimates for the power waveforms shown in 
Fig. 2.3 for cases 1 and 2, square wave pulses will again be assumed as 
shown in Fig. 2.b (d, e). 

Basically, the calculation of the failure levels will be done in 
three parts: 

1. Calculation of the temperature distribution throughout 
the semiconductor due to the first pulse. 

2. Calculation of the temperature distribution at the end of 
the cooling period. 
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SINGLE PULSE (Actual data) 



TIME (Microseconds) 









SINGLE PULSE (Actual data) 



TIME (Microseconds) 






TIME (Microseconds) 









3. Calculation of the temperature at the junction due to 
the second pulse. 

A. First Pulse Heating 

The equation for the temperature divided by the power density 
is given by 


T_ 

P 

A 


2 



erfc 


2/.ct/(pC p ) 


( 11 ) 


A plot of the relative temperature versus distance from the junction for 
various pulse durations is shown in Fig. 2.10. The distance at which the 
temperature is down by l/e from the junction temperature is marked. 


B. Temperature After Cooling 

A good approximation to the cooling can be found by using the 

equation 

T = {erf — - -- + erf - - * X -- — ■ } . (12) 

2/Kt/(pC) 2/Kt/(pC ) 

P P 

The above equation is valid for the case where tho region -a < x < a 
is initially at constant temperature T q and the region |x| > a is 
initially at zero. 

From the temperature distribution of Part A (Fig. 2.10), we 
will estimate a value for 'a' for various pulse durations and calculate 
the normalized temperature distribution for various cooling durations. 
Table II gives the values of 'a' chosen for various heating-pulse 
durations. 
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Figure 2.10 







u 


TABLE II 

Pulse Length 


a 

(Microseconds) 


(Microns) 

0.01 


o.u 

0.10 


1.0 

1.00 


h.O 

10.00 


10.0 

100.00 


hOO.O 


Figures 2.11, 2.12, and 2.13 are plots of the normalized temperature 
versus distance from the junction for various cooling times. The plots 
are for initialized heating due to 0.1, 1.0, and 10.0 microsecond pulses 
respectively. 

C. Second Pulse Heating 

Finally now we need to calculate the temperature rise due to 
the second pulse. Because th temperature of the device is no longer a 
constant throughout the device, we can no longer use the equation of 
Part A to calculate the temperature. In order to do this exactly, one 
needs to go back to the original heating equation and derive a new 
expression using a continuous plane source at th'e junction and include 
the actual temperature distribution as found in Part B. However, a 
reasonable estimate of the temperature due to the second pulse may be 
found by using equation H with a time duration equal to the actual pulse 
duration plus an "effective time" of the first pulse and cooling period. 

This "effective time" will be equal to the time required to heat the 
material to the temperature after the cooling period due to a power level 
of the first and second pulses. First, we need to derive the expression 
for ' fc e ff ec ^i ve ) 2nd find equation for the failure level due to 

both pulses. We will use f fi to denote the cooling factor at the junction 
as found from Figs. 11 to 13. 
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TEMPERATURE OF SEMICONDUCTOR REGION FOR VARIOUS COOLING TIMES 



DISTANCE FROM JUNCTION 
Figure 2.11 









TEMPERATURE OF SEMICONDUCTOR REGION FOR VARIOUS COOLING TIMES 












TEMPERATURE OF SEMICONDUCTOR REGION FOR VARIOUS COOLING TIMES 









Figure 2.14 


Using equation h, we can write the equations 

for the various 

temperatures ps shown in Fig. lh as 


[T -T J = 7 K t 1/2 , 

1 o A ’ 

(13) 

IT JS] ml K t 1/2 

1 2 V A K eff * 

(HO 

and 

f'f _«n ] = £ k t 1/2 

U 3 *o J A 3 * 

(15) 

where 


t 3 = t + t eff and K s (tucpC ^)" 1 ^ 2 

(16) 

Also 


T 0 ■ f o 'W + T o 

(17) 

Equation lU then becomes 


{f [T.-T ] + T -T } = 7 K t^ 2 

(18) 






Solving for t 


u 


f c [T i- T o ]2 


(19) 


From equation 13 

Then (19) becomes 


P 1/2 + 

1 A o 


f 2 if- K T 1/2 + T -T ] 2 
+ _ c L A _ o & _ -2 _ 

eff “ ,P ^2 " c 


( 20 ) 


( 21 ) 


•o-f' 


. )W2 = E, 

eff ^ A 


( 22 ) 


p /fticpCp [T^ -Tj I' 


- 1/2 


(23) 


From equation 23, we can now calculate our failure levels for 
a double pulse with a cooling time between pulses. Note that for no 
cooling, f Q = 1 and (l + f 2 ) 1 / 2 = /2, which makes equation 23 reduce 
exactly to equation 7, which we had derived previously. An example of how 
the cooling period affects the failure level is shown for the transistor 
case in Fig. 15 with the calculated values given in Table III. The 0.1 
and 10 microsecond points differ slightly relative to the 1 microsecond 
point because the length of the heated region was chosen at the distance 
where the temperature was down by 36 . 8 $ (i.e., 1 ) for the 1.0 microsecond 
pulse and at the distance where the temperature was down by U6.8$ for the 
0.1 and 10.0 microsecond points. 
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TABLE III 


Double 

Pulse Failure 

Condition 

Value of K of Eq. 10 

1.0 microsecond pulses 

No cooling 

219 

1.0 microsecond pulses 

1.0 microsecond cooling 

277 

1.0 microsecond pulses 

2.0 microsecond cooling 

291 

0.1 microsecond pulses 

0.1 microsecond cooling 

907 

10.0 microsecond pulses 
10.0 microsecond cooling 

90.7 


DISCUSSION 

Although we feel that the multiple pulse analysis discussed is 
generally valia, there may be some devices whose failure levels may 
differ from those predicted. 

For large area diodes we have noted that pulses of long duration 
and large powers can be dissipated without damage until a certain power 
level - time duration is reached. If the power level is then raised 
slightly the device fails in very short times. We believe that this 
phenomena is due to the diode initially conducting over most of the 
junction area until a certain temperature and electric field gradient is 
reached. Then a small hot spot or breakdown occurs in a very localized 
region. When this occurs, a heavy current is channeled through the hot 
spot and the device quickly fails. What happens then when pulse trains as 
shown in Fig. 2.2 are applied to such a device? In this case, a steady 
power level and unidirectional voltage is not maintained. The lower power 
level and the change in voltage polarity as the voltage wave goes into 





L. 

the second pulse may sufficiently change the conditions for breakdown 
and hot spot formation that the device might not fail as readily. 

On the other hand, a very rapid change in applied voltage polarity may- 
result in some very high local field gradients in the junction and break¬ 
down conditions might be enhanced. However, in either case, drastic 
deviations from the estimated failure levels are not expected. 

Another precaution one should note is on cases 1 and 2, where we 
have assumed that reverse conduction does not occur. For small signals 
this is certainly true; however, for the power levels we are discussing, 
the first conduction pulse will cause very high junction temperatures. 

At such temperatures, the reverse voltage breakdown characteristic is 
drastically changed. Conduction may then actually occur during the 
reverse cycle. The momentary change in breakdo'vn voltage cannot be 
noted from our single pulse experiments since the device has sufficient 
time to cool before a measurement can be made on the transistor curve 
tracer. 

^ C ONCLUSIONS 

It is felt that the foregoing analysis shows that one does not expect 
drastic changes in the failure levels due to a series of several pulses. 
The general order of magnitude failure level still holds. Even for the 
worst case of assuming continuous power flow for the sum duration of all 
the pulses, the results are still of the same order of magnitude. Any 
cooling shifts the curves back toward the single pulse failure level. 

Also, as shown in Fig. 9, the estimated failure level for a triple pulse 
is actually still within ohe spread of experimental data points for a 
single pulse. 

The worst case condition described by equation 7 is probably adequate 
to describe all multiple pulse conditions if one is only interested in an 
order of magnitude estimate of the failure level. 
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SECTION 3. 


INTRODUCTION 

I. Semi-Empirical Threshold Power Failure Equatio n 

A final re-evaluation of test data for all the Pershing semiconductor 
damage tests produced the results shown in Tables I and II of this section. 
The significant numbers in these tables are listed under column heading 
^(watt-sec 1 ^ 2 ). These numbers allow semi-empirical threshold power failure 
levels to be obtained from the relation 


p _ a .-i/2 

P thres " K * 


Threshold power failure levels will be obtained in watts when pulse widths, 
t, are expressed in seconds. This is the semi-empirical equation (based on 
the thermal failure model) which permits order of magnitude estimates of the 
failure level by adjusting the Kr constant to best fit the data. Presented 
in this manner, the semiconductor damage results are in the proper form for 
circuit-response-to-transient analyses and further interpretation of the data 
is unnecessary. 

II. Threshold Power Failure Levels for Devices When Damage is Achieved 
The general experimental procedure to determine the threshold power 
failure level of semiconductor junctions consists of voltage pulsing the 
junction until failure occurs. "Failure" is arbitrarily defined to be a 
15% reduction in transistor DC current gain (Hpg) or a similiar degradation 
in diode zener voltage. From this data, plots of power versus pulse duration 
are obtained. Based on these plottings, a semi-empirical threshold power 
failure curve is drawn. This curve reflects the reverse failure points, no 
failure points, and the uncertainty in the data for each point. The failure 
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curve is described by the equation 


P = K 1 t- 1/2 . 

From the test results obtained, it can be predicted with a high 
degree of confidence that power levels an order of magnitude higher than 
those levels estimated by the threshold power failure equation will produce 
device failure; conversely, power levels an order of magnitude lower than 
those indicated by the equation will not cause device degradation. 

III. Threshold Power Failure Levels for Devices When Damage is Not Achieved 

Device failure for a number of large-area diodes and transistors was 
not possible within the limits of the BDM pulser. For these devices, 
threshold power damage levels were estimated from junction areas and the 
"standard" P/A versus time curves shown in Figure 3.1. These curves are 
defined by P/A = K^/A t -1 ^, where the constant, K^/A, represents the 
average value obtained from diodes or transistors which were damaged. 

The resulting K /A constants are 

~ = ( 0 . 551 * x 10 3 ) watt-sec^^/cm^ for diodes 

1 *' 1/2 

= (0J17 x 10 3 ) - for transistors 

cm 

Multiplying the appropriate average K"*VA value above by the junction 
area of an undamaged device yields an "estimated" K. Having this estimated 
value, a threshold power failure curve can be drawn for the diode or 
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transistor not damaged. 

A slight departure from the prediction of the threshold power failure 
equation was denoted for the large area devices. As a result, -it is 
anticipated that a higher power failure level is estimated by this equation 
for those undamaged, large area devices. This undetermined error is especially 
true for components having junction areas measuring in tenths of square centi¬ 
meters . 


IV. Device Damage Plots 


Plots of power versus pulse duration for each device tested are contained 


in this section. The solid line on the plot represents the threshold power 
failure curve; the dotted line drawn on the plot represents the estimated 
threshold power failure curve defined in Paragraph III. 
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SECTION k. 


NOHSEMICONDUCTOR DEVICE TESTING 


INTRODUCTION 

Various electrical components (relays, capacitors, coils, etc.) vere 
evaluated for damage vulnerability due to high voltage transients- In 
general, two areas were investigated for these nonsemiconductors: (l) 
device degradation resulting from the applied pulses (dielectric break¬ 
down, changes in electrical responses or characteristics, etc.); and (2) 
internal arcing occurring during the pulse application. Virtually 
all of these device types proved to be "hard" for the maximum levels 
attainable with the BDM pulser. 

Table I of this section summarizes the test results of the 
nonsemiconductor testing. This tabulation includes all the Pershing 
devices tested. 

TEST PROCEDURE 

Each device tested was multiple pulsed (10 times) across their terminals. 

In general, the pulser output was increased in 3 steps to its maximum output 
voltage (UCO vo.lts, 700 volts, and 1000 volts) for each tested component. The 
applied pulse width was 8 p-seconds. Voltage and current photographs were 
taken for each of the selected pulse amplitudes and from these recordings the 
test data was reduced. 

Degradation to any device was established by comparing various measurements 
taken by an impedance bridge before and after an applied pulse series. Using 
the bridge, the basic electrical characteristics could be measured for the 
various samples tested. Additionally, seme components, such as motors and 
relays, had their rated potentials applied to them after the damage tests 
and observed for abnormal behavior during actual operation. 





VACUUM TUBE DAMAGE TESTING 


A cursory damage study was made to establish the vulnerability of 
vacuum tubes to high voltage transients. The complete experimental damage 
test results for tubes 6BCU, 6BX7, and 5876 is presented in Table II of this 
section. Although a brief study was conducted, the results did indicate that 
permanent damage to vacuum tube devices is possible with high voltage pulses. 



TABLE I. NONSEMICONDUCTOR DEVICES (PERSHING) 


Device Description 

Pershing Number 

Failure 

Devices Tested 

Brake Clutch 

1061*0ll*3 

0 

2 

Capacitor 

5910 - 060-1188 

0 

1 

Capacitor 

5910-669-3625 

0 

1 

Capacitor 

5910-8ll»-61»U9 

0 

1 

Capacitor 

5910-823-1068 

0 

1 

Capacitor 

KF223KM* 

0 

1 

Coil 

595C-679-9539 

0 

2 

Connector 

10509676 

0 

1 

Filter 

MIS 17279/1-1 

0 

1 

Filter 

MIS 171*9/6-1 

0 

1 

Filter 

1101*0553-1 

0 

1 

Inductor 

5950-81*5-6927 

0 

1 

Motor 

MIS 17103/1-3 

0 

1 

Potentiometer 

106071*95-7 

0 

1 

Potentiometer 

106085 }* 8-1 

0 

1 

Potentiometer 

1062071*95-11 

0 

1 

Potentiometer 

1101*0271-3 

0 

1 

Relay 

10607305 

0 

1 

Relay 

10630138 

0 

1 

Relay 

1063011*2 

0 

1 

Relay 

1101 * 001 * 7-1 

0 

1 

Relay 

1101*061*7 

0° 

2 

Switch 

1063009-1 

0 

1 

Switch 

1063009-3 

0 

1 

Transformer 

MIS 17088/2-1 

0 

1 

Transformer 

51,1,-9715-002 

0 

i 

Transformer 

5950-856-6285 

0 

1 • 

Transformer 

1101*01*99 

0 

r. 

Transformer 

1101*051* 1»-1 

0 

1 

Transformer 

1101*051*5-1 

0 

1 

Transformer 

1101*0590-1 

0 

1 


* Manufacturer's Number (WES CAP) 

0 Internal Arcing During Applied Pulses. Intermittent arcing at approximately 
500 volts, constant arcover at approximately 700 volts. 






TABLE II. VACUUM TUBE DAMAGE TEST RESULTS 


TUBE 

6BX7 1 
5876 2 
6bc4 3 


Notes: 


NUMBER 

TESTED 

FAILURE 

GRID-CATKODE 

VOLTAGE 

PULSE 

WIDTH 

3 

0 

+ 1KV 

8us 

2 

0 

+ 1KV 

8ps 

3 

2* 

+ 1KV 

3us 


* Transconductance (gm) decreased approximately 35# for Loth 
damaged devices. 

Medium - mu twin triode. 

2 

Pencil type UHF high-mu triode. 

3 Medium - mu triode. 





APPENDIX A 


JUNCTION MODELING 

THEORETICAL MODEL F OR JUNCTION FAILURE DUE TO TEMPERATURE 

When a voltage pulse is applied to a semiconductor device in the 
reverse direction, the principal voltage drop is across the junction. 

An approximate thermal model for this situation is a plane source at 
the junction x = 0 in an infinite medium whose thermal conductivity and 
specific heat are functions of temperature. For silicon, the thermal 
conductivity is a sensitive function of temperature, e.g., from 1.56 
watts/cm-°K at 300°K fco 0.310 watts/cm-°K at 1000°K [2l»]; and since 
large temperature changes are anticipated, the change in thermal 
conductivity must be considered. 

When the thermal conductivity and specific heat are functions of 
temperature, the general one-dimensional heat equation is [ 25 ]: 


8_ 

8x' 


dl \ _ dl 

n 


0 


(1) 


where k = thermal conductivity [watt/cm-°K]; p = density [grams/cm^]; 

0^ = specific heat [joules/gm-°K]; t * time [sec]; T = temperature [°K]; 
and x = distance [cm]. 

By a change of variables, 


9= /^ aT ’ 

O 

equation can be written as [23]. 


( 2 ) 


8j0 1_ 86 

9x 2 - 8(6) 8t 


0 


(3) 


This equation is similar to the familiar heat equation except that a is a 
variable which depends on 6. For the semiconductor failure problem, actual 
junction geometries vary from device to device, and the exact junction 
mechanisms are not necessarily consistant. Because of these uncertainties. 
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certain simplifying assumptions are justified. Therefore, temperature 
variations for the semiconductor failure model will "be described by 


u 


3 2 T 1 3T 

3x 2 "“ at 


(10 


where a is a constant. Corrections for possible over-simplifying assumptions 
can be made by using an ’effective’ or time-weighted average of k over the 
temperature range of interest. The weighting will be done by considering 
the time interval in each temperature interval considering the change in 
1/2 

temperature is given by AT = At 7 (which we will find is the form of our 
solution). 

The solution of equation (4) for an infinite medium ,/ith an instantaneous 
plane source of strength Q, at t = 0, parallel to the plane x = 0 and passing 
through x' with zero initial temperature is given by 


Q -(x-x') 2 /4at, 

(W ) 1/2 


(5) 


where the heat liberated per unit area is QoC^. For a rectangular heat source, 
of duration t, Q = Q(t') = q = constant, and the temperature rise is 


T 


-U~x ') 2 

t l»a(t-t’) 
(4ira) 1/2 / (t-t’) 1/2 


dt' . 


( 6 ) 


Performing the integration, 


T = q (—)^ 2 e Ucrt" - erfc/1 


kxiL\ 

\2fsr j 


(7) 


Making the substitution x = 0 ana x' = 0, then equation (7) fits the 
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requirements of the thermal model for semiconductor heating at the junction, 
and 


G 


T(0) = q (r-)' 


t_vl/2 


( 8 ) 


represents the junction temperature increase caused by a pulse of width t 
which heats the junction at a rate q. The heat per unit area per unit time 
is given by the power per unit area divided by the specific heat and the 
density, i.e.. 


<1 


P 1_ 

A P C p ‘ 


(9) 


If the initial junction temperature is T^, then T -*■ AT = [T-T^], and equation 
(8) becomes 


[T-T.] 


P 

A 




( 10 ) 


Experimental guidelines are obtained from equation (10) if temperature 
extremes, within which junction failure is deemed likely to occur, can be 
estimated. Defining T m to be the junction failure temperature, the power 
per unit area as a function of time is given by 

| -lf55c p [T^-Tj] t' 1/S (11) 

Three cases are of interest. The constants for these cases are: 

p = 2.33 gm/cm 3 , C p = 0.7566 joules/gm-°K, < S ± = 25°C, and k 
( effective) = 0.526 watts/cm-°K (Case I), and k (effective) = 0.306 watts/cm-°K 
(Cases II and III). (For time in microseconds, P/A is given in kilowatts/cm 2 .) 

Case I - Heating from room temperature 25°C (~300°K) to a 

temperature of 675°C (9^8°K). Equation (11) is then: 

P/A = 1109. t" 1/2 


{ 
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Case II - Heating from room tempo', -«.<,ure 25°C (~300°:v) to the 
melting temperature of .ilicon Al5°C (i688 £ 'K), 

(excluding phase char;.*.- Equation (ll) is thent 

P/A = 1809.Y t' U ' : 

Case III - Same as Case II except that the current is assumed to 

pass through local hot spots and the power is dissipated 
in one-tenth of the total area. The failure temperature 
is the hot spot temperature. 

Semiconductor failure levels should fall within these limits, and plots for these 
three cises should give an indication of the reliability of experimental data. 

As previously mentioned an effective value of thermal conductivity was 
used in the calculations above. The thermal conductivity of silicon is neither 
a constant nor a linear junction of temperature as indicated in Fig. A-l. 

Values of the thermal conductivity k for 100 C K intervals are given in Table A-l. 

It should be noted that the first and last values are given for 50°K intervals. 
Since AT = C^'^Tt", setting AT = 100 for a full interval for a unit time will 
give C 1 = 100. Setting AT = 50 fo? the first step will then given t = lA. 

The till'e differential At^ repress iL.- a relative time increment between intervals 
weighted with respect to the time during which power is delivered to the junction, 
whereas t.. represents the relative time with reference to the total pulse duration. 
The value of k^ is the thermal conductivity at the midpoint of the intervals and 
represents an average value in the temperature range indicated. 

The table values given can be used to determine the effective thermal 
conductivity in an arbitrary temperature range according to the formula: 


eff E A t. 


( 12 ) 


EXAMPLE 1. Determine the effective thermal conductivity for silicon 
in the temperature range from room temperature (300°K) to a terminal 
temperature of 950°K. 










Table A-l Values of Thermal Conductivity for Silicon 


t(°k) 

AT. 

1 

t. 

3. 


i 


300-350 

2 J 1 

l/*4 

1/4 

1.360 

0.34 

350-450 

In 

2 C 1 

9/1 1 

2 

1.050 

2.10 

450-550 

5. Q 

2 C 1 

25/4 

4 

0.800 

3.20 

550-650 

K 

49/4 

6 

0.640 

3.84 

650-750 

K 

81/4 

8 

0.520 

4.16 

750-850 

11 c 

2 c i 

121/4 

10 

0.430 

4.20 

850-950 

13 c 

2 C 1 

169/4 

12 

0.356 

4.27 

950-1050 

15 „ 

r u i 

225/4 

14 

0.310 

4.34 

1050-1150 

lie 

2 C 1 

289/4 

16 

0.280 

4.48 

1150-1250 

12 . c 

2 c i 

361/4 

18 

0.261 

4.70 

1250-1350 

21 c 

2 c i 

441/4 

20 

0.248 

4.96 

1350-1*450 

53 c 

2 

529/4 

22 

0.237 

5.21 

1*450-1550 

2f c 

2 C 1 

625/4 

24 

0.227 

5.45 

1550-1650 

21 c 

2 C 1 

729/4 

26 

0.219 

5.69 

1650-1700 

28 c 

2 C 1 

784/4 

13-3/4 

0.216 

2.97 
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SOLUTION: 


Sum of column At.^ from 300°K to 950°K = 22.21 
Sum of column At., from 300°X to 950°iC = k2.2 c - 


P. 2.21 
K eff " ¥2.25 


0.52b y/cm-°K. 


EXAMPLE 2. Determine the effective thermal conductivity of silicon in 
the temperature range from room temperature (300°K) to l,7'00 o K (melting 
point of silicon). 

SOLUTION: 

Sum of column At^ic. from 300°K to 1, 700°K = 60.01 
Sum of column At^ from 300°K to 1,700°K = 196 


K - 60-01 
eff " 196 


0.306 w/cm-°K. 


The thermal conductivity versus temperature for germanium is given in 
Fig. A-2. A table of values for germanium similar to that discussed for 
silicon is given in Table A-2. The effective value for the thermal conductivity 
of germanium for the temperature interval 300°K to 1210°K (melting point of 
germanium) was calculated to be 0.202 watts/cm-°K. Solving equation 11 
for germanium using the constants p = 5.33 gm/cm , ” 0.3093 ,1oules/gm- 0 K, 

= 25°C, = 937°C, and k effective = 0.202 watts/cm-°K, the power per 

unit area for failure is given by 





u 


Table A~2 Values of Thermal Conductivity for Germanium 


T(°K) 

AT 

t 

300 
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Figure A-2. Thermal Conductivity of'Germanium 




Theoretical Voltage and Current Calculations 
1. Reverse Polarity Calculations 

The calculation of the "breakdown voltage in the reverse direction 
has been adequately described for certain conditions- In the reverse 
direction, two different microscopic mechanisms - occur: 

1. Avalanche breakdown and 

2- Zener breakdown 

The normal VI characteristic for a p-n diode, I = I [e^^"-l] gives 
the correct current in the reverse direction [V negative] and then I = -I q 
when the absolute magnitude of -V is much much greater than kT/q until some 
maximum reverse voltage is reached. Then the equation no longer applies 
and the reverse current increases very rapidly. The physical phenomenon 
which cause-, the reverse breakdown is called the avalanche effect. It is 
a carrier multiplication effect which occurs in the p-n junction when a 
sufficiently large voltage is applied in the reverse direction. It is a 
non-destructive mechanism in itself. The reverse current characteristic is 
microscopically described as follows. As the holes [or electrons] which 
constitute the stauration current I q move through the junction region, they 
gain energy from the electric field and then give it up when they collide 
with the lattice. For smaller electric field intensities, all the energy 
is given up to the lattice with each collision with the lattice and the hole 
continues through the junction with no increase in current. As the electric 
field intensity increases, the hole gains more energy per mean free path 
length and due to quantum-mechanical considerations, does not give up all 
of its energy with each collision. Hence, the hole keeps gaining energy 
with each mean free path until it has sufficient energy to ionize the 
lattice upon collision, thereby, releasing more free hole-electron pairs. 
This multiplicative action continues and the avalanche current is produced. 



If the current is limited by circuit resistance or by time durations, 
there is no damage to the device. The avalanche breakdown is the normal 
breakdown mechanism in junctions with light doping and wide junction regions. 

The Zener breakdown is the field-induced direct transition of carriers 
from the valence band to the conduction band. This breakdown (essentially 
field emission) occurs when the electric field strength reaches values of 
approximately 10^ volts/cm cr higher.' As the doping level increases, the 
mean free path length decreases so with narrow junction widths, a particle 
cannot gain sufficient energy to cause ionizing collisions, but the field 
is strong enough to cause direct transition. Thus, again the current 
sharply increases from its saturation level due to the Zener breakdown 
which again is non-destructive if the current is limited. (This same 
effect is the one that changes the shape of the forward current characteristic 
and gives a negative resistance region for tunnel diodes.) 

Again, there has been considerable theoretical work done in calculating 
the reverse breakdown voltage. These calculations agree well with experimental 
evidence if the physical parameters of the junction are known. 

Although reverse breakdown voltages can be theoretically calculated 
it does not appear to be worthwhile to do this for the study of failure 
levels. In order to make accurate calculations the physical parameters 
of the junction need to be accurately known. The needed information is 
often not readily available from the manufacturer. The reverse breakdown 
voltage for the device can be easily measured experimentally. As long as 
the reverse current is limited during the test there is no damage to the 
device. The breakdown voltage determined experimentally is more accurate 
and gives an exact breakdown level for the particular device to be tested 
later for failure level. 

2. Forward Polarity Calculations 

The calculation of current and junction voltages for forward conduction 
has been quite well established for the lower current levels. The current 
calculated from I = Io(e^^^-l] matches closely with experimental values. For 
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large currents these relationships are no longer valid. At large currents 
there is an appreciable IR voltage drop across the bulk resistance of the 
diode. Kano and Reich [26] have discussed a model for calculation of transient 
parameters when the forward current levels are quite high. However, their 
model is based on a step impulse of current rather than voltage. Their results 
indicate that the equation above is still valid for quite large currents provided 
that the voltage value used in the equation is the actual voltage across the 
junction rather than the terminal voltage. Thus the current is given by 


I = lo 


wheiTG y — \j _ y 

junction terminal bulk resistance. 

A plot of the V-I characteristics are as shown in Figure A-l. 



Figure A-l 

In order to make temperature rise calculations, the voltage and the 
current across the junction must be known. 
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Since the current rises as an exponential function of the junction 
voltage, this voltage can be considered to be a constant in calculating 
the power dissipated in the junction. Attention is now given to the 
temperature calculations in order to determine the failure levels. The 
power can be calculated from P = V.Io e ^ V j^ kT for high current levels. 
To perform a more accurate failure calculation, the heat dissipated in 
the bulk resistance must be considered since it can be appreciable at 
high current levels. 
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APPENDIX 3 

TRANSISTOR BASE-BUTTER JUNCTION PROTECTION 
AGAINST HIGH VOLTAGE TRANSIENTS 3Y THE USE OF DIODES 

I. GENERAL 

Previous work performed by the laboratory indicates that the transistor 
junction most susceptable to high voltage transients is the base-emitter. 
This is primarily because of its low zener voltage and small cross-sectional 
area. 

This report describes several means of protecting this junction 
against damage with the use of high peak inverse voltage diodes and zener 
diodes. The advantages and disadvantages of these protection devices will 
also be discussed. 

Since the laboratory has obtained a good collection of data on the 
2N2222 transistor, this device was used in the experimental test conducted 
for this study. 

H. 2N2222 CHARACTERISTICS AND DAMAGE LEVELS 

A. Characteristics 

The base-emitter junction of a typical 2N2222 is rated at a 
reverse bias breakdown voltage ( v BE q) equal to 6 volts DC. The mean value 
of the devices tested was 7.8 volts DC with a spread being from 7.0 to 
8.8 volts DC. The emitter cut-off current which flows from the emitter to 
the base is rated at a maximum value of 10 nanoamperes DC with a 3 volt 
DC reverse bias present on the base emitter junction. Beyond the zener 
knee, the dynamic resistance of the junction is 7.5 ohms about 200 
milliamperes. 

In the forward direction, the junction exhibits a resistance of 
2.3 ohms at a forward current of 0.5 amps DC and I .56 ohms at 1 ampere. 

The dynamic resistance in both cases is approximately 0.8 ohms. 

B. Damage Levels -. 

Typical failure points for the 2N2222 in the reverse and forward 
direction are shown in Figure 1. Each point represents a to'cal loss of the 
base-emitter diode characteristic. No attempt is made here to indicate 
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degradation in zener potential or reverse leakage current. The points 
indicated are representative of the actual peak junction dissipation in the 
reverse direction of a single high voltage transient applied to the junction. 
Junction current and voltage data was obtained with a fast trace recording 
camera. 

III. PROTECTION AGAINST REVERSE TRANSIENTS 

A. Series Diode 

The base-emitter junction is capable of conducting large amounts 
of peak power because of its very low zener potential (6-10 volts). The 
zener potential can be increased without degrading the normal operation 
of a circuit by inserting a high peak inverse voltage diode in series with 
the emitter of the transistor (see Figure 2A). This combination will not 
draw reverse current until the zener potential of the additional diode is 
exceeded. 

A test'was run with a 1R645 diode (which has a rated PIV equal to 
275 volts DC) in series with the emitter of a 2N2222. Where rhe 2N2222 
usually conducts 25 amperes with an applied 1KV pulse, the current drawn 
for the combination was less - chan 80 milliamperes. The junction was not 
damaged by conduction times as long as one microsecond. 

The advantage of this technique is that high peak inverse voltage 
diodes have very low capacitance, high back resistance and a relatively 
low forward drop so as not to degrade the circuit's normal performance. 

The major disadvantage, of course, is that the diode does not 
offer any forward transient protection. 

B. Shunt Zener Diode 

Another method of reducing the dissipation of the base-emitter 
junction of the transistor is to shunt the current around the junction 
and dissipate the power in a stronger device (see Figure 2B). 

The zener diode will protect the junction in the reverse transient 
case by conducting in its normal forward direction. Care must be used to be 
sure that a temperature compensated zener is not used as this would require 
a very high voltage to overcome the PIV of the compensating diode. 





A. SERIES DIODE {REVERSE PROTECTION) 



B. SHUNT DIODE (REVERSE PROTECTION) 



A+ 


C. SHUNT DIODE WITH SERIES Rg (FORWARD PROTECTION) 



Figure 2. Protection Circuits Against Transients. 
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For the forward transient case, the base-emitter junction would 
conduct before the zener even with a low voltage zener. High voltage 
transients would cause the zener diode to avalanche but not until the 
base-emitter junction had already heavily conducted large amounts of current, 
damaging the junction. 

Tests conducted show that indeed the shunt zener did protect 
against reverse damage transients but failed to protect against the forward 
transients. 

Another disadvantage of the shunt diode protection circuit is the 
large amount of capacitance that the diode, presents to the driving source. 

A typical value of 1 volt reverse bias capacitance for a low voltage 
zener is on the order of 500 picofarads. This additional capacitance, 
depending upon the applications, may be highly degrading upon the operation 
of the circuit. 

IV. PROTECTION AGAINST FORWARD TRANSIENTS 

A. Shunt Zener Diode with Series Base Resistance 

The reason the shunt zener diode failed to protect the base- 
emitter junction against forward transients is because of the zener impe¬ 
dance being higher than the base-emitter saturation resistance. Increasing 
the saturation resistance of the base-emitter junction by inserting a small 
resistor in series with the base will cause the zener to handle the 
majority of the current. The small resistor will not effectively degrade 
the normal performance of the circuit (see Figure 2C). 

Data obtained from experimental tests of this configuration 
indicated that with a 12 ohm base series resistance the 3.3 volt zener diode 
conducted approximately 32 amperes while the base-emitter junction con¬ 
ducted only 0.5 amperes with a 1 kilovolt pulse applied to the configuration 
in the forward direction. No damage wa^* done to either of the semiconductor 
deivces. 

Removal of the zener diode and testing with only the series base 
resistance indicated the expected failure. 
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The shunt sener diode in this configuration has the same dis¬ 
advantages mentioned previously in Section III-B. The large shunt capaci¬ 
tance may he detrimental to the normal operation of high-speed circuitry. 
V. OTHER PROTECTIVE DIODES INVESTIGATED 

Thyristors, specifically four-layer trigger diodes and switching 
tunnel diodes were briefly investigated as possible protective devices 
because of their high switching speed. A review of the possible applica¬ 
tion of these devices indicated that the limited dissipation capabilities 
and self susceptability of these devices to transients of this type 
prohibi: their use. The limited number of devices which have higher power 
dissipation capabilities do not appear to be any harder. 





















APPENDIX C' 


THRESHOLD POWER FAILURE CURVES FOR DEVICES 2112222, MC26JG, MC355G, ALT) 
SNT311. 

Concurrent with the Pershing semiconductor damage study, additional 
devices were tested as part of BDM "in-house" study program. Extensive 
damage testing of the 2112222 silicon KPN transistor was conducted and 
the results from this device were used to assist in the validation of 
the thermal failure mode}. Also, three integrated circuits, MC267G, 
MC355G, and SN73U, were studied. The damage information from the three 
I. C.'s served as a comparison to that data gotten from the discrete 
transistors and diodes tested for the overall program. 

The following pages of this appendix contain the threshold power 
failure curves for the in-house study devices. Note that for the 2N.2222 
tz*ansistor, curves are also presented for conditions where the case 
temperature was elevated before pulse application to the B-E junction. 
The complete results of all these tests are included in m ables I and II 
of Section 3. 
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